Nonequilibrium Precursor Model for the Onset of Percolation in a Two-Phase System 
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Using a Boltzmann-likc equation, we investigate the 
nonequilibrium dynamics of nonperturbative fluctuations 
within the context of Ginzburg-Landau models. As an illus- 
tration, we examine how a two-phase system initially prepared 
in a homogeneous, low-temperature phase becomes populated 
by precursors of the opposite phase as the temperature is in- 
creased. We compute the critical value of the order parameter 
for the onset of percolation, which signals the breakdown of 
the conventional dilute gas approximation. 
PACS number(s): 05.70.Ln, 05.70.Fh, 64.60.Ak 

The process by which an ordered low-temperature sys- 
tem approaches disorder as its temperature is increased 
has far-ranging applications in many areas of physics. 
In systems which allow for low-temperature symme- 
try breaking, such processes would describe symmetry 
restoration. There is considerable overlap between treat- 
ments found in the condensed matter literature within 
the context of Ginzburg-Landau models p],p| and those 
found in the high-energy physics literature through the 
use of temperature-corrected effective potentials |3fl, al- 
though clearly there are several crucial differences as well 
j|. In this letter, we would like to explore an issue which 
is of interest to both areas, namely, how to describe the 
dynamics of nonperturbative thermal fluctuations in sim- 
ple systems, modeled by the Ginzburg-Landau model. 

It is well-known that fcrromagnets will become para- 
magnetic above a certain critical temperature. It is also 
well-known that such emergence of disorder is due to 
the nucleation of ferromagnetic domains of the oppo- 
site phase ||. The dynamics of the domain interfaces, 
as well as their growth, is of interest in many diverse 
areas, from materials science to particle physics to cos- 
mology, even if some systems require more complicated 
order parameters. Examples include the recent exper- 
iments on Bose-Einstein condensation in dilute atomic 
gases H and the study of the growth of the condensate 
0; ultra-relativistic heavy-ion collision experiments 
formation of topological defects both in the laboratory, 
via pressure quench experiments , or in the early Uni- 
verse ]f0[; a nd the nematic-isotropic transition in liquid 
crystals ]I1|]. Diverse as these systems are, they all have 
one feature in common; their change in behavior is due to 
the onset of nonequilibrium conditions, which are poorly 
understood. 

In studying these kinds of problems one usually starts 



with mean-field theory, or some of its microscopic ver- 
sions, such as the equilibrium one-loop approximation in 
field theory or the Hartree approximation. In these ap- 
proximations, localized, high amplitude fluctuations are 
neglected and replaced by an average interaction of the 
system with the thermal environment Q. However, this 
approximation breaks down as the system approaches 
criticality and these fluctuations become more important. 

It has long been recognized [Q that in order to fully 
understand the dynamics of a given system one must go 
beyond the mean field approximation. One approach is to 
invoke time-dependent renormalization group techniques 
0. Here, we would like to propose an alternative ap- 
proach, based on the dynamics of large-amplitude fluc- 
tuations, from which we can examine the nonequilibrium 
properties of the system. As we will show, our approach 
is valid up to the onset of percolation, which is known to 
occur before criticality for 3-dimcnsional systems. 

Let us start by considering a standard Ginzburg- 
Landau model where local fluctuations about homoge- 
neous equilibrium have the free energy 



F(<j>,T) 



|v^| 2 + y(^) 



(i) 



with V((j>) = V + a(6 - I)0 2 /2 + A0 4 /4, where a, b, V 
and A are (positive) constants and 9 is the temperature 
ratio T/T c . For convenience, we have added the constant 
term Vo to fix the minima of the free energy density at 
T < T c at zero, (0) = ±y/a(l - 6)/\ = <j>±, which then 
gives that Vq = a 2 (I — 6>) 2 /(4A). The constants a and 
b can easily be scaled away and F(<f>, T) can be made 
dependent only on the temperature ratio 6 and on the 
coupling constant A. 

We choose to study the dynamics of the fluctuations 
as the system is heated from its T = state, where it 
is in one of its ordered states, say 0_, to a temperature 
T < T c , focusing on its evolution to a final equilibrium 
state determined by a time-independent value of the or- 
der parameter (</>) at temperature T. We model the fluc- 
tuations away from the initial equilibrium state as Gaus- 
sian shaped, spherically-symmetric configurations with a 
core value tf>c an d radius R. These precursors are also 
called subcritical bubbles and treatments involving them 
have been successfully used in many other contexts be- 
fore pa] . By expressing the amplitude of the fluctuation 
as 4>a = 4>c ~ 4>- 1 t ne fluctuations are parameterized as 
4>f(r) — (f>A exp(— r 2 /R 2 ) + </>_. [The Gaussian satisfies 



f 



the physical boundary conditions - regularity at the ori- 
gin and asymptotic approach to the background, while 
costing very little free energy. See Ref. |l3| for details.] 

Since we are interested in fluctuations which can probe 
the other available free-energy minimum, their ampli- 
tudes can be easily determined by the condition that 
4>f(r) represents unstable fluctuations inside the (— )- 
vacuum phase. We then simply have, from the sym- 
metric double well potential used in Eq. ([!]) (note that 
small-amplitude fluctuations are already incorporated in 
the mean- field approximation to the free-energy density) , 
that those fluctuations with (pc > </>max = are the ones 
probing the (+) phase. As for R, we allow for fluctua- 
tions larger than the correlation length, i? m i n = £(T), 
where f(T) = [V"(</>±)] -1/2 , consistent with the nat- 
ural coarse-graining scale dictated by the continuous 
free energy. Substituting <f>f(r) in ([l]) we obtain the 
free energy barrier for a fluctuation with amplitude 4>a 
and radius R as F f (R,(f> A ,T) = (3V2tt 3 ^<P 2 a /8)R + 

^/ 2 4> 2 A [V2V"{4>-)I% " + 

y/W"'{4>-)4)A/^ + V ( - 4 \<p-)(t) 2 A /l92}R 3 . We will refer 
to the fluctuations with 4>c > </>max as "(+)-phase fluc- 
tuations" and the background phase as the "(— )-phase". 

We next study the dynamical evolution of these fluc- 
tuations. For this we use the Boltzmann-like equa- 
tion derived in Ref. Jl4| for the distribution function of 
(+)-phase fluctuations of radius R and amplitude (j>A, 
f+(R, 4>a, t), which satisfies the equation 



df± 
dt 



V] ~dR +( 1 -7)G-^+-7G+_ , (2) 



where the first term in the rhs incorporates the col- 
lapse of subcritical domains, which we approximate as 
having constant velocity v = dR/dt. In a forthcom- 
ing publication we will show that this is a valid ap- 
proximation for most of the interesting range of bub- 
ble radii. The second term describes the thermal nu- 
cleation of fluctuations of the (+)-phase inside the (— )- 

phase, with nucleation rate G >+, while the last term 

describes the inverse process, with rate G_| For a de- 
generate potential these two rates are the same, which 
we express in terms of the free energy of the fluctu- 
ations, Ff, as given by a standard Gibbs distribution: 
G(<f>A,R) = G-^+ = G+^_ = AT 4 exp(-F//T), where 
A is a constant. Note also that from Eq. (g), detailed 
balance imposes that the ratio be constant, which 
will be taken as a free parameter in our model; it can be 
determined for specific models, as shown in Ref. ]T^ |. In 
fact, the ratio must depend on dynamical quanti- 
ties which are, in principle, expressable in terms of the 
only two parameters in the model free energy, A and 8, 
which must control heat diffusion and fluctuations dy- 
namics and can be mapped to specific applications. 7 in 
Eq. (||) is the total fraction of volume in the (+)-phase, 
defined by Q 



+00 ,+oc 4 r3 

d<j>c dR — - f+(R,<i>A,t) . (3) 



Note that from our initial condition at t = 0, we have 
j(t — 0) = and, in the asymptotic equilibrium regime 
at temperature < T < T c , < 7 oq < 1/2. Eq. (§), 
together with (j^) , is an integro-differential equation that 
we numerically solve for 7. The result is shown in Fig. 1 
for different values of temperature and A = 1 . 
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FIG. 1. The volume fraction 7 as a function of time for 
A = 1 and A/\v\ = 1. 

7 C q can be computed from Eq. (|^) by setting the time 
derivative term in the lhs to zero. The resulting expres- 
sion for 7 oq , using Eq. (||), is 7 oq = 1/(1 + 21), where 
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Writing j(t) = j eq B(t), and using Eqs. (||) and (||), we 
find an equation for B(t), after integrating all terms in 
Eq. (|) by L + °° d^ c dR (AttR 3 /3) , 



^W + iL B ( t )-iL = o 

dt 7oq 7oq 



(5) 



where T, the total volume-integrated nucleation rate, is 
given by 

r+00 ,+00 4 ^3 

T = / d^ c dR^—G(cbc:R) • (6) 

The differential equation (^) has a simple solution, given 
by B(t) = 1 — exp(— t/r), where r = 7 cq /r is the equi- 
libration time-scale. Therefore, the analytical solution 
for 7(i) is j(t) = 7 oq (l — e~*/ r ). This solution fits very 
well the numerical solution for 7 shown in Fig. 1. In 
Fig. 2 we compare the theoretical and numerical results 
for the equilibration time-scale r, as a function of the 
temperature, for two different values of A; the agreement 
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is quite striking. The results in Fig. 2 reveal a peak in 
the equilibration time-scale. For small temperatures, the 
equilibration time grows continuously until it reaches a 
maximum at the temperature T max (A). This is reminis- 
cent of the critical slowing down behavior characteristic 
of critical phenomena, although we cannot recover the 
discontinuity at the critical point with our simple model. 
What we do provide is a dynamical picture of the ap- 
proach to criticality, which we expand on below. As the 
temperature increases, a larger fraction of the volume of 
the initial state in the (— )-phase is populated by fluctu- 
ations of the (+)-phase; also, as the free-energy barrier 
decreases with increasing temperature, these fluctuations 
will become larger in size. Hence, their equilibration 
time-scales grow with increasing T, as displayed in Fig. 
2. For temperatures larger than T max , the equilibration 
time decreases, vanishing at the mean-field critical tem- 
perature T c . This is due to the fact that at T > T max the 
system is described by a free-energy density with a single 
minimum at = 0; thus, the true critical temperature is 
not at T c . This is in accordance with what is expected 
when corrections beyond the mean field to the potential 
are taken into account Jl5| and explicitly seen in large- 
scale Langevin simulations performed on the lattice [^6| . 
From Fig. 2 our model predicts the values T max ~ 0.79 T c 
for A = 1 and T max ~ 0.97 T c for A = 0.1. These results 
are dependent on the ratio which entails the mi- 

croscopic details of a given model under study. Physical 
lower bounds on this ratio will be discussed below. 
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FIG. 2. The equilibration time r for A = 1 and A = 0.1, 
with A/\v\ = 1. The dots are the numerical results from Eq. 
([]) and the lines the theoretical result, T t h C or = 7cq/r. 

As the temperature is increased and the volume- 
density of (+)-phase fluctuations grows, the system will 
eventually reach an instability point beyond which do- 
mains of the (+)-phase will grow by percolating with 
their nearest neighbors [jl7|. The question is at what 
temperature such percolating instability occurs. In order 
to answer this question, we take full advantage of our 
dynamical model. [A preliminary approach can be found 



in Ref. H].] Since correlation- volume fluctuations have 
the smallest free-energy barrier, they will be statistically 
dominant. In order to model the percolation- instability, 
consider a domain of the (+)-phase of correlation-length 
radius £. There are three main processes that can change 
its volume: shrinking due to its surface tension; growth 
due to the thermal nucleation of another (+)-phase do- 
main of radius R just outside it, which accounts for 
a change of volume AV = 47r[(£ + 2R) 3 - £ 3 ]/3; and 
thermal destruction of the correlation size fluctuation 
due to inverse nucleation, that changes the volume by 
AV' — 47r£ 3 /3. We thus arrive at an approximate equa- 
tion describing the rate of change of Vt : 

-±^-\ v \^e+ d^ c dR^-G(4> c ,R)AV 

-/ dcjjc-^G^c^AV'^^e^s. (7) 

In Fig. 3 we show the numerical solution for the ef- 
fective velocity, w G ff, as a function of temperature. The 
temperature for which v e g > 0, T pcrc (A), indicates the 
onset of percolation for correlation- volume fluctuations. 
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FIG. 3. The effective velocity of correlation length fluctua- 
tions, Weft = (4-7r£ 2 ) dVf/di as a function of temperature for 
A = 1 and A = 0.1. A/\v\ is again set to one. 

Our numerical results give the values T porc ~ 0.72 T c for 
A = 1 and T pcrc ~ 0.96 T c for A = 0.1. T pcrc obtained 
with equation (Q) is very close to T max given before. As 
percolation sets in, Eq. (0) begins to underestimate the 
growth of fluctuations and the further development of the 
system. T pcrc then determines the limit of validity of Eq. 
(0), or the breakdown of the dilute gas approximation, 
beyond which coalescence of phase fluctuations begins to 
be of importance. Nevertheless, Eqs. (Q) and (^) describe 
quite well the dynamics until the onset of continuous per- 
colation as well as the equilibrium properties of the sys- 
tem. (Note that continuous percolation, as opposed to 
lattice percolation, is very poorly understood, and only 
within simple 2-dimensional mathematical models, such 
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as the Boolean- Poisson model Jl7|.) 

Finally, it is important to test the validity of this model 
with respect to the calculation of the free energy of the 
fluctuations. It is clear that as the free energy Ff for fluc- 
tuations drops below fc^T we no longer can distinguish 
them from simple thermal noise, which then becomes sta- 
tistically dominant; the description of the nucleation of 
large-amplitude fluctuations with rate G becomes mean- 
ingless. Using the temperatures T max and T pe rc, we can 
set a lower bound on the value of This is shown in 

Fig. 4 for the case of A = 1 , where we show how T max and 
Tpcrc change with ^4/|v|. The condition Ff/T\x n = 1 ap- 
plied to the fluctuations of lowest free energy, <f) — 4>max 
and R = f , gives T n > 0.87 T c for A = 1 and T„ > 0.99 T c 
for A = 0.1. 
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FIG. 4. # max = T max /T c and per c = T pclc /T c as a function 
of v4/H, for A = 1. The dotted line is denned by the condition 
Ff /T\t„ = 1, from which we obtain the bound A/\v\ > 1CF 2 . 

Summarizing, we have presented a simple model based 
on the dynamics of phase fluctuations that is able to pro- 
vide a dynamical description of how a continuous ordered 
system described by a Ginzburg-Landau free energy ap- 
proaches its percolation instability. The model allows 
us to compute the temperature for the onset of perco- 
lation, which signals the breakdown of the conventional 
dilute gas approximation, offering also a novel way to 
estimate the actual critical temperature in a Ginzburg- 
Landau system, that can be read from the results of Fig. 
2. Furthermore, the model studied here, despite its sim- 
plicity, exhibits a dynamical picture of symmetry restora- 
tion and the breakdown of mean-field theory observed 
both numerically and analytically, without recourse to 
large-scale numerical simulations. The model can eas- 
ily be extended to different systems including inhomoge- 
neous nucleation, or systems outside the Ising universal- 
ity class. We expect to report soon on these applications. 
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